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ABSTRACT 

We present the Spitzer Space Telescope Infrared Spectrograph spectrum of the Orion A protostar 
HOPS-68. The mid-infrared spectrum reveals crystalline substructure at 11.1, 16.1, 18.8, 23.6, 27.9, 
and 33.6 /mi superimposed on the broad 9.7 and 18 /im amorphous silicate features; the substructure 
is well matched by the presence of the olivine end-member forstcrite (Mg2Si04). Crystalline silicates 
are often observed as infrared emission features around the circumstellar disks of Herbig Ae/Be stars 
and T Tauri stars. However, this is the first unambiguous detection of crystalline silicate absorption 
in a cold, infalling, protostellar envelope. We estimate the crystalline mass fraction along the line- 
of-sight by first assuming that the crystalline silicates are located in a cold absorbing screen and 
secondly by utilizing radiative transfer models. The resulting crystalline mass fractions of 0.14 and 
0.17, respectively, are significantly greater than the upper limit found in the interstellar medium 
(< 0.02-0.05). We propose that the amorphous silicates were annealed within the hot inner disk 
and/or envelope regions and subsequently transported outward into the envelope by entrainment in 
a protostellar outflow. 

Subject headings: circumstellar matter — infrared: stars — stars: formation — stars: individual 
(FIR-2) — stars: protostars 



1. INTRODUCTION 

Since their discovery in the envelopes around oxygen- 
rich evolved stars, crystalline silicates (e.g., pyrox- 
ene [(Mg, Fe)Si0 3 ] and olivine [(Mg, Fe) 2 Si0 4 ]) have 
been detected in the circumstellar environments of both 
pre- and post-main-sequence stars, comets , and ul- 
traluminous infrared galaxies (jHennin and ref- 

erences therein). However, silicates in the interstel- 
lar medium are alm ost entirely amorphous in structure 
([Kemper et al.l l2005). Crystalline silicates are thought to 
be produced by thermal annealing or evaporation and re- 
condensation in warm environments; experimental stud- 
ies indicate tha t such processes require temperat ures of 
- 800-1200 K (|Harker fc Deschl 120021 [C^l200l . 

In pre-main-sequence stars, crystalline silicates are 
often detected in the warm inner disk region (< 
1 AU), where the temperature is sufficient to anneal 
amorphous silicates (iBouwman et aT1 [2008: Sar gent et al.l 
l2009t lOlofsson et al.l 120101 and references therein). In 
contrast, there are only a few claims of crystalline sili- 
cates in protostars, which are still surrounded by cold, 
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infallin g envelopes a nd th us in an earlier evolutionary 
stage. ICiardi et al.l (|2005l) attributed short-wavelength 
emission features in the protostellar binary SVS 20 to 
crystalline silicates. However, such emission must origi- 
nate in the warm inner regions near the central protostar 
and may be produced in a circumstellar disk as in pre- 
main-sequence stars. Conversely, crystalline silicates in 
the cold, infalling envelope would result in absorption 
features. Shallow absorption near 11.3 has been de- 
tected toward several lo w-mass protostars (jBoogert et al.l 
12004 iRiaz et al.l 12009ft . The associ ation of t his fe ature 
with crystalline silicates is debated: iLi et al.1 (|2007l ) sug- 
gest instead that the inclusion of a H2O ice mantle on 
silicate grains resu lts in a weak sho ulder near 11 /tm. 
In support of this, IRiaz et ail (|2009l ) find a strong cor- 
relation between the 11.3 /zm feature and the H2O ice 
column density. 

In t his Letter, we pre sent the Spitzer Space Tele- 
scove (iWerner et "all 120041 ) Infrared Spectrograph (IRS; 
iHouck et al.1 12004T ) spectrum of the dee ply embedded 
proto star HOPS-68 (also known as FLR-2: [Mezger et al.l 
1990?). Th e source is located in the Orion Molec- 
ular Cloud 2 regi o n (O MC-2; ICatlev et al.l Il97l 
iPeterson &: Megeathl [2008, and references therein), 
northward of t he Orion Nebula, at an adopted distance 
of 414 ± 7 pc l|Menten et al.l 120071 ). We report the first 
detection of a complex of absorption features which can 
unambiguously be attributed to crystalline silicates in 
the cold envelope of a protostar. 

2. OBSERVATIONS AND DATA REDUCTION 

The 5-36 /im spectrum was obtained toward the 
Spitzer 4.5 /im position of HOPS-68 (a = 05 h 35 m 24!3, 
5 = -05°08'30" [J2000]) on 2007 March 27 (Spitzer AOR 
20838656) with the short- (SL; 5.2-14.0 /mi) and long- 
wavelength (LL; 14.0—36.1 /im) low-resolution (A/AA = 
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Figure 1. Top panel: Spitzer-IRS 5—36 fira spectrum of the protostar HOPS-68, with uncertainties (gray lines) and spectral features 
indicated. The feature at 9.66 fim is a background-subtraction artifact, arising from the H2 S(3) emission line near the source. Bottom 
panel: The spectrum of the subtracted background. Although PAH features and H2 lines are detected, crystalline silicate features are not 
apparent. 



60—120) IRS modules at each of the two nominal nod 
positions, one-third of the way from the slit ends. The 
spectrum was extracted from the Spitzer Science Cen- 
ter (SSC) S18.7 pipeline basic calibrated data using the 
Advance d Optimal Extracti on in the SMART software 
package (Hig don et al.l 12004]) follow ing the methods de- 
scribed in lLebouteiller et al.l(|2010D . The LL spectra and 
first order SL spectra were background-subtracted by fit- 
ting a degree zero polynomial to either side of the emis- 
sion profiles. For the SL module, the background emis- 
sion was scaled up (down) by 30% before taking the nod 
difference for the first (second) nod positions. 

The resulting nods were then averaged to obtain the 
final spectrum, and the spectral uncertainties are esti- 
mated to be half the difference between the two inde- 
pendent spectra from each nod position. We estimate 
the spectrophotometric accuracy of the final result to be 
5%. 

Additional 70 and 160 /im imaging was performed 
on 2010 September 28 (observation IDs 1342205228 
and 1342205229) with Herschel' s Photodetector Ar- 
ray Camera and S pectrometer (jPilbratt et al.l 120101 : 
IPoglitsctTe t al. 201(3) as part of the Herschel Orion Pro- 
tostar Survey (HOPS). An 8' field was repeatedly ob- 
served eight times using two orthogonal scanning direc- 
tions and a scan speed of 20" s _1 , giving a total obser- 
vation time of 3290 s. The data were proc essed with 
"Method 1" described in lFischer et~aH (|201Q[ ). 



Aperture photometry was obtained using a 9.6" and 
12.8" aperture at 70 and 160 /<m, respectively, with 
the background annulus extending from 9.6—19.2" and 
12.8—25.6", respectively. An aperture correction was 
determined from the encircled energy fraction provided 
by the PACS consortium (private communication). The 
photometric accuracy is dominated by calibration un- 
certainties, which are 10% and 20% at 70 and 160 xxm, 
respectively. 

3. ANALYSIS 

The Spitzer-IRS 5—36 /xm spectrum exhibits copi- 
ous silicate and ice absorption features (Figure [T]) that 
are typically found in the circumstellar environments of 
embe dded protostars (|Boogert et al.l 2004; Watso n~et al.l 
120041 ). The main ice constituents are H2O (6.0 and 12.0 
H, CH3OH (6.8 and 9.9 /xm), CH 4 (7.7 /xm), NH 3 (9.0 
lira), and CO2 (15.1 /xm). In addition, the spectrum re- 
veals two gas-phase organic molecules: C2H2 (13.7 /xm) 
and HCN (14.0 /xm); these have been previo usly observed 
toward massive protostars (|An e t al. 2009]). The promi- 
nent silicate absorption features at 9.7 and 18 /xm are 
due to the Si-0 stretching and O-Si-0 bending vibra- 
tional modes, respectively. 

The spectrum also reveals a complex of narrow features 
near 11, 16, 19, 24, 28, and 34 /im. These features are 
typic ally detected in emission around the disks of T Tauri 
stars (Sar gent et al.ll2009l : lOlofsson et al.| [2010) and have 
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Figure 2. The isolation of crystalline silicate features in the 
spectrum of HOPS-68. Top panel: The Spitzer-YRS spectrum 
with the adopted second-order polynomial unabsorbed continuum 
(gray line). Second panel: The derived optical depth spectrum 
compared to the amo rphous silicate profile of GCS-3 (gray line; 
Chiar & Ticlcns 2006). Third panel: The residual optical depth 
spectrum, after subtracting the scaled amorphous silicate profile of 
GCS-3, compared to a laborat ory spectrum of pure crystalline H2O 
ice at T = 140 K (gray line; Hudgins et al. 1993). Bottom panel: 
The residual optical depth spectrum after subtracting the scaled 
H2O ice profile. Arbitrarily scaled mass absorption coefficients for 
forsterite (solid gray line) and enstatite (dotted line) are shown 
for comparison; thes e we re calculated using optical constants from 
Sogawa et al. (2006) and Jager ct al. (1998), respectively. 



been attributed to the presence of crystalline silicates. 
Moreover, this complex of features has previously only 
been detected in absorption to ward heavily obscu red ul- 
traluminous infrared galaxies (jSpoon et al.ll2006h . 

3.1. Cold Absorbing Screen Method 

To investigate the presence of crystalline silicate ab- 
sorption, we approximate the envelope as a cold ab- 
sorbing screen again st a bright infrared continuum (e.g., 
iBoogert et al.ll2008T ) . The continuum is approximated as 
a second-order polynomial that is fit to the regions of the 
spectrum that appear devoid of absorption features: 5.5, 
7.5, and 31.6 fim. 

The optical depth spectrum, r(A), is derived assum- 
ing F£ hs = i^oixt e -T(A) 5 where F obs and F cont are the 

observed and continuum flux densities, respectively. To 
account for the amorphous silicate contribution, the op- 
tical depth spectrum is compared to the silicate ab- 
sorption profile of the Galactic Center source GCS-3 
(jChiar fc Tielen s 2006) , which is commonly used to rep- 
resent "standard" interstellar silicate absorption. The 
polynomial continuum fit was offset to a greater than 
observed flux at 7.5 and 31.6 /zm; the applied offsets en- 
sured an observed optical depth greater than or equal to 
the combined GCS-3 and H2O ice optical depth for all 
wavelengths. 

The top panel of Figure [2] shows the result after ap- 
plying the offsets to the polynomial continuum fit. The 
GCS-3 optical depth was then constrained to fit the ob- 
served optical depth spectrum at 9.3—9.5 /jm by mini- 
mizing 

x = ? MAO ' (1) 

where T a (Xi) and er (Ai) are the observed optical depth 
and its corresponding uncertainty at the ith wavelength, 
respectively, r m (Ai) is the tabulated optical depth at the 
ith wavelength for GCS-3, and a is the scaling parameter 
of the fit. The scaled amorphous component is shown in 
the second panel of Figure [2] Next, the residual optical 
depth spectrum is compared to the labor atory spectrum 
of pu re crystalline H 2 ice at T = 140 K (Hudgins ct al. 
[1991 . which was calculated for spherical grains in the 
Rayleigh-limit (A. C. A. Boogert, private communica- 
tion). The H 2 ice profile is fit to match the residual op- 
tical depth spectrum between 12.6—13.0 fan (third panel 
of Figure [5]). The crystalline H2O ice optical depth spec- 
trum provides a good fit to the depth of the broad H2O 
libration mode, and thus is preferred over the laboratory 
spectrum of pure amorphous H2O ice. 

Subtraction of the scaled H2O ice profile from the op- 
tical depth spectrum isolates features near 11, 16, 19, 
24, 28, and 34 /jm, as illustrated in the bottom panel of 
Figure [21 The peak positions, widths, and strengths of 
the residual features are inconsistent with the strongest 
features of the py roxene end- member enstatite (MgSiOa ; 
IJager et al.l 119981 ). Conversely, the peak positions and 
widths of the residual features at 11.1, 16.1, 18.8, 23.6, 
27.9, and 33.6 /im coincide with the strongest of those 
in the mas s absorption coeffici ents (MACs) for forsterite 
(IVfeSiO^ ISogawa et al.ll2006ft . The relative strengths of 
the 18.8 and 23.6 ^m features appear weaker than those 
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Table 1 

Crystalline Mass Fraction 



A 

H 


^cr (-V)) a 

[10 3 cm i g" 1 ] 


Tex (^p) b 


N a /(N CI +N &m ) c 
[Hudgins] 


N CI /(N CI +N aui r 
[Warren] 


11.1 


6.0 


0.90 (0.12) 


0.15 (0.02) 


0.13 (0.02) 


16.1 


1.5 


0.26 (o.oi) 


0.16 (o.oi) 


0.20 (o.oi) 


18.8 


3.8 


0.22 (0.02) 


0.06 (o.oi) 


0.09 (o.oi) 


23.6 


4.2 


0.23 (<o.oi) 


0.06 (<o.oi) 


0.07 «o.oi) 


27.9 


1.2 


0.15 (<o.oi) 


0.13 (<o.oi) 


0.11 (<o.oi) 


33.6 


2.4 


0.43 (0.02) 


0.17 (o.oi) 


0.15 (o.oi) 



Note. — Uncertainties based on statistical errors in the Spitzer-IKS spectrum 
only. Peak optical depth of the 9.7 fim amorphous silicate feature measured from 
the GCS-3 fit: r am (A p ) = 3.05. 
a Peak mass absorption coefficients for forsterite. 
b Peak optical depth of forsterite. 

c Cr ystalline mass fraction s deri ved from Equa tion 2 using H2O optical constants 
from Hudgins ct al. (1993) and Warren (1984), respectively. 



calculated for forsterite; a similar mismatch between 
the observed and c alculated strengths was reported by 
ISpoon et al.l J2006). Nonetheless, we conclude that the 
features observed at 11.1, 16.1, 18.8, 23.6, 27.9, and 33.6 
/xm are caused by the presence of the mineral forsterite. 

The crystalline mass fractions are inferred from the 
peak optical depths of the 9.7 /xm amorphous, T am (A p ), 
and the 11.1, 16.1, 18.8, 23.6, 27.9, and 33.6 (im crys- 
talline, r cr (A p ), silicate absorption features using 



Tcr(Ap)Kam(A p ) 



A a 



r(Ap) K am(Ap) + T am (A p )K C r(Ap) 



(2) 



where N cr and N aia are the dust mass column densities 
for crystalline and amorphous silicates, respectively. The 
peak MACs for forsterite, K cr (A p ), are tabulated in Ta- 
ble [I] and we r e calc ulated using optical constants from 
iSogawa et all ([2006D f 33~2l . Adopting «; aro (A p ) = 3.5 x 
10 3 cm 2 g -1 as the peak MAC for amorphous pyroxene 
([Sargent et al.l [20061 ). we find crystalline mass fractions 
ranging from 0.06—0.17, with a median value of 0.14. 
These crystalline mass fractions are significantly greater 
than the u pper limit found in the interstellar medium (< 
0.02-0.05; lKemper et al JI2005I: ETet al Jl2007n . A factor- 
of-three variation in the crystalline mass fraction may 
be due to uncertainties in the MACs for forsterite, un- 
certainties in the subtracted optical depth spectra for 
the amorphous silicate and H2O ice contributions, and 
finally the assumption that all the silicate material is 
in the cold absorbing screen. Furthermore, assuming 
more realistic grain shapes for forsterite may reduce the 
discrepancies among the crystalline mass fractions, but 
we assume a continuous distribution of ellipsoids (CDE; 
iBohren fc Huffman 1983) for simplicity. 

The analysis is repeated u sing the H2O ice optical 
constants from iWarrenl ([1984D . The results show that 



the mass fraction of forsterite does not depend on the 
adopted H 2 ice optical constants (Tabled]). 

3.2. Radiative Transfer Method 

To account for the presence of both warm emitting 
material in the inner envelope and cold absorbing ma- 
terial in the outer envelope, we now analyze the spec- 
trum using a radiative transfer model. The model is 
bas ed on the method in troduced in lKenvon et al.l (119931) 
and lCalvet et al.l (|1994D . and later refined in lOsorio et al.l 



( 2003) and iFurlan et ail (|2008l) to include the sheet- 
collap se models developed by IHartmann et al.l (1994, 
1996). The sheet-collapse model uses the solution for the 
collapse of a self-gravitating, isothermal, infinite sheet 
initially in hydrostatic equilibrium; this model results in 
a relatively flat spectral energy distribution (SED) over 
the mid-infrared wavelength region. 

A central protostar and disk are responsible for heating 
the surrounding envelope. The disk is modeled as a flat, 
irradiated, steady accretion disk that is optically thick at 
all wavelengths. The most important input parameters 
include: the total system luminosity (L = £ s tar+-E'acc)> 
the fraction of the total luminosity that arises from the 
star (77 s tar), the centrifugal radius (R c ), the inner (i? m i n ) 
and outer (i? max ) envelope radii, the inclination angle to 
the line-of-sight (i), the reference density at 1 AU (pi), 
and the degree of asphericity (77 = R ma ^/H, where H is 
the sc ale height of the original sheet; see IHartmann et"aLl 

mil). 

We assume a si milar envel ope dust mod el as de scribed 
in IQsorio et all ([200l and IFurlan et al] ([200l . The 
grain size distribution follows the standard power-law 
distribution n(a)da = a _3,5 c?a, with 0.005 /im < a < 
0.3 /xm. The fractional abundance with respect to the 
mass of gas for t he ice and dust sp ecies is similar to 
those proposed by IQsorio et al.l ([20031 ): £sii = 0.003, C gra 
= 0.0025, Ctro = 0.000768, Cice = 0.0012, and Cco 2 = 
0.00039 for amorphous silicate, graphite, troilite, H 2 
ice, and CO2 ice grains, respectively. Details of the 
calculat ion and references for optical properties c an be 
found inlD'Alessio et ail (|2001| ). [Osorio et al.l([2003D . and 
IFurlan et al.l 



In addition, we include a CDE grain shape distribution 
of forsterite in the Rayleigh-li mit. We adopt optical con - 
stants for (Mg, Fe) 2 Si0 4 from Huffman fc_Stapp| (|1973ft . 
Mgi.q FeruSiCU from IZeidler et al.l (I2011D. Mg9SiC>A 
from ISogawa et alJ ([2006D . and Mgi gFeo iSiC>4 from 
IFabian et all ([200ll) for the 0.1-0.16, 0.24-2, 4-100, 
and 100—8000 /mi wa velength regimes, res pectively. The 
optical constants from lFabian et all (120011) were modified 
to match the values from lSogawa et alJ (|2006[ ): the real 
part n was subtracted by a constant, while the imagi- 
nary part k was multiplied by a scalar. The MACs for 
forsterite were computed from the complex indices of re- 
fraction for each crystalline axis, and the three resultant 
values were then averaged together to account for ran- 
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Figure 3. Left panel: Best-fit model compared to the observed SED of HOPS-68. The composite SED includes the Spitzer-IRS spectrum, 
PANIC i<" s -band flux (J. J. Tobin, private communication), IRAC 3.6, 4.5, 5.8, and 8.0 /im and MIPS 24 (im fluxes (Meg eath et al., in prepa- 
ratio n), Herschel-PACS 70 and 160 fim fluxes (F„ = 6.9 and 23.3 Jy, respectively), and SCUBA 450 and 850 (im fluxes ( Johnstone fc Ballvl 
1999). The SCUBA estimates represent peak fluxes assuming a point-like source. The model {blue line) includes contributions from the 
central protostar {magenta line), the disk {green line), and the envelope {orange line); the displayed disk and central protostar components 
are extinguished by the envelope. In addition, a gray-body {dash-dotted line) is also included to simulate external heating of the outer 
envelope. Right panel: Close-up of the 5—40 pia region with IRS uncertainties {gray lines) and forsterite wavelength positions indicated 
{dotted lines). Envelope models with increasing crystallinity {violet lines; 0.0, 0.13, and 0.20, respectively) are shown for comparison 
against the best-fit model {blue line; 0.17). The mismatch near 8—9 pra between the adopted "astronomical" silicates and the observed 
amorphous silicate feature has been noted in the literature fHcnning 20 1(3). 

(1.3 Lq) protostar with a flattened (r) — 2.0), relatively 
high density envelope (pi = 5.7 x 10 -14 g cm~ 3 ) and a 
small centrifugal radius (R c — 0.5 AU), and is viewed at 
an inclination of 41°. (Assuming a central mass of M* = 
0.5 Mq, the reference density suggests a mass infall rate 
of M = 7.6 x 10" 6 M Q vr- 1 : iKenvon et all [i"99l . A 
small centrifugal radius may be the result of a very slowly 
rotating molecular cloud prior to collapse or may indicate 
a very short elapsed time since the onset of collapse. 

The best-fit SED model has a forsterite mass frac- 
tional abundance of Cfor = 0.00063 relative to gas and 
adequately fits the depth of the 11.1, 16.1, 18.8, 23.6, 
27.9, and 33.6 /im forsterite features (right panel of 
Figure 131 ■ The best-fit mass fractional abundances for 
crystalline and amorphous silicates (Cfor/(Cfor + Csil) = 
0.00063/0.00363) imply a crystalline mass fraction of 0.17 
This result is consistent with the cold absorbing 
screen method; the slightly larger mass fractions from 
the radiative transfer model may result from forsterite 
emission features generated in the warm (> 100 K) in- 
ner envelope, which must then be absorbed by the cooler 
outer envelope. 

4. DISCUSSION 

We have identified forsterite (Mg 2 Si04) absorption in 
the cold, infalling envelope of a protostar, and estimate 
a crystalline mass fraction of 0.14 and 0.17 from the 
cold absorbing screen and radiative transfer methods, re- 
spectively. These values are three to four times greater 
than the upp er limit found in the interstellar m edium 
(< 0.02-0.05: Kemper et all 120051: ILi et all 120071) . The 
good agreement between the crystalline mass fractions, 
determined using two different methods with different 
optical constants, demonstrates the robustness of this 
result. Even if we adopt the discrepant mass fraction 



Table 2 

Best-Fit Model Parameters 

Parameter Value 
Sheet-Collapse 
L [Lq] 1.3 

??star 0.3 

n 2.0 

R c [AU] 0.5 

Rmin [AU] 0.39 

i?max [AU] 7000 

pi [lO" 14 g cm" 3 ] 5.7 

i [dcg] 41 

9 [deg] a 18 

CforACfor + Csil) ■• 0.17 

Gray-Body 

I., K 2L5 

rioo 0.012 

P 1.17 

xl d o_*_ 

a Cavity semi-opening angle. 

domly oriented grains (Boh ren fc Huffman! Il98l . 

The SED and best-fit model for HOPS-68 are shown 
in the left panel of Figure [3] Following the procedures of 
iFurlan et al.l ()2008l ) . the model parameters were adjusted 
by visual examination to yield a good fit over the entire 
SED with more emphasis given to the mid-infrared. To 
account for the observed increase of flux beyond 70 /im, 
we include a s imple gray-body function (dash-dotted line; 
IGordonlll995[ ) to simulate the effects of external heating 
on the outer envelope; the dust temperature (Td), emis- 
sivity index (/3), and optical depth at 100 fim (ti 00 ) are 
determined using a nonlinear least-squares fitting rou- 
tine. The best-fit parameters are tabulated in Table [2] 

The SED modeling indicates a moderately luminous 
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of 0.06 estimated from the 18.8 and 23.6 /xm forsterite 
features (Table [TJ, the observed crystalline mass fraction 
still exceeds the upper limit for the interstellar medium. 

Crystalline silicates are o bserved in emission toward 
the disks of young stars fe.g-. ISargent et aLll2009D . where 
the temperature in the inner disk region is adequate 
for annealing; however, it is unlikely that the observed 
forsterite absorption arises in the disk. While absorption 
features are produced in the atmospheres of rapidly ac- 
creting disks, such observation s typically show fe atures 
longward of 10 /im in emission (|Green et al.1120061 ) . Fur- 
thermore, since the contribution of the envelope is more 
than a factor-of-ten greater than the attenuated disk 
component (Figure [3]) , an absorption feature in the disk 
would require an unphysical depth to produce an ob- 
served depth of 20% in the combined SED. 

If the envelope is primarily material transported from 
the interstellar medium through gravitational collapse, 
what process is responsible for producing the observed 
crystalline grains? For a slowly rotating, flattened en- 
velope, much of the infalling material is deposited close 
to the central protostar, where it is then readily heated 
to temperatures sufficient for annealing. But, because 
forsterite absorption is detected in the 20—30 /mi regime, 
the absorbing grains must be located in a relatively low 
temperature (< 100 K) region (r > 15 AU). It is highly 
unlikely that the amorphous silicates are annealed in situ 
by irradiation, since raising the dust temperatures at dis- 
tances > 15 AU above the glass t ransition temperature 
(T glass = 990 K; iSpeck et alJl200l requires luminosities 
> 10 4 Lq. Therefore, an alternative scenario for the 
production of crystalline silicates in the outer envelope 
is needed. 

We propose that the amorphous silicates were annealed 
within the hot (> 1000 K) inner disk and/or envelope 
regions (< 0.4 AU) and subsequently transported out- 
ward into the surrounding cold envelope material by en- 
trainment in a protostellar outflow. (An outflow associ- 
ated w ith the protostar was detected by Williams et al.l 
(2003).) Alternatively, the amorphous silicates in the 
envelope may have been annealed in situ by outflow- 
driven shocks; however, it is not clear that such shocks 
could heat grains to 1000 K without de stroying them via 
sputtering or gain-grain collisions (e.g., iDraine fc McKed 
1993; iNeufeld fc Hollenbach|[l99l . 

If these mechanisms result in short annealing times in 
which chemical equilibrium is not achieved, they would 
produc e the obser ved dominance of forsterite over en- 
statite l)Gailll200"4T ). While further analysis is needed to 
better understand the origin of the observed forsterite, 
its presence suggests that crystalline silicates detected in 
the outer disks of young stars may have been deposited 
during the protostellar infall phase. 
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